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O Abstract 

^ In the outer envelope of the Sun and in other stars, difFeren- 
O tial rotation and meridional circulation are maintained via 

the redistribution of momentum and energy by convective 

piotions. In order to properly capture such processes in a 

numerical model, the correct spherical geometry is essential. 

Jn this paper I review recent insights into the maintenance of 
' ^ 'mean flows in the solar interior obtained from high-resolution 

^simulations of solar convection in rotating spherical shells. 
I The Coriolis force induces a Reynolds stress which transports 
2 angular momentum equatorward and also yields latitudinal 

variations in the convective heat flux. Meridional circula- 
te f;ions induced by baroclinicity and rotational shear further 

redistribute angular momentum and alter the mean stratifl- 
I cation. This gives rise to a complex nonlinear interplay be- 
^ tween turbulent convection, differential rotation, meridional 

circulation, and the mean speciflc entropy profllc. I will de- 

scribe how this drama plays out in our simulations as well as 

in solar and stellar convection zones. 
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1 Introduction 

Axisymmetric flows are a key component of virtually all so- 
lar dynamo models. Differential rotation is thought to be the 
principal mechanism by which toroidal field is continually re- 
generated from poloidal field. In the language of mean-field 
dynamo theory, this is known as the fi-effect and the Sun may 
be classified as an a-fl dynamo (e.g. Ossendrijver 2003; Char- 
bonneau 2005). Furthermore, many recent dynamo models, 
known as flux-transport models, suggest that global axisym- 
metric circulations in the meridional plane (latitude-radius) 
may account for the observed migration of active regions to- 
ward the equator during the course of the solar activity cycle 
and may thus determine the cycle period (e.g. Wang et al. 
1991; Choudhuri et al. 1995; Dikpati & Charbonneau 1999; 
Dikpati & Oilman 2006). 

The internal rotation profile of the Sun is well established 
from helioseismology (Thompson et al. 2003) and as we have 
seen at this meeting, there is ample evidence for the presence 
of differential rotation in other stars (see other papers in these 
proceedings as well as previous work by Barnes et al. 2005 
and Reiners 2006). In the solar convection zone the angular 
velocity decreases monotonically by about 30% from equator 
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to pole with nearly radial contours at mid latitudes. Regions 
of stronger radial shear exist near the top and bottom of the 
convective zone, the latter known as the solar tachocline. 

The meridional circulation is more challenging to detect 
but photospheric measurements and local helioseismic inver- 
sions have revealed a systematic but variable poleward fiow 
of 15-20 m s~^ in the upper solar convection zone (e.g. Hath- 
away 1996; Haber et al. 2002; Gonzales-Hernandez et al. 
2006). 

Differential rotation and meridional circulation in the so- 
lar envelope and in other stars are maintained by convection. 
The highly turbulent nature of solar and stellar convection 
makes numerical modeling a difficult challenge but contin- 
inuing advances in supercomputing technology have enabled 
increasingly realistic 3D simulations. In this paper I'll review 
what our most recent convection simulations reveal about the 
nature of mean fiows in the solar envelope and how they are 
maintained. 

2 Turbulent Solar Convection 

High-resolution simulations of solar convection exhibit intri- 
cate flow structures and persistent mean flows as illustrated 
in Figure 1. The simulation shown was carried out using 
the ASH (Anelastic Spherical Harmonic) computer code de- 
scribed by Clune et al. (1999) and Brun, Miesch & Toomre 
(2004). The ASH code solves the three-dimensional equations 
of hydrodynamics (or magneto-hydrodynamics) in a rotating 
spherical-shell geometry under the anelastic approximation. 
The simulation shown in Figure 1 extends from the base of 
the convective envelope at 0.71i?, where R is the solar ra- 
dius, to 0.98i?, within about 14 Mm of the photosphere. Be- 
yond 0.98i?, granulation and supergranulation occur in the 
Sun which cannot be resolved in a global model (yet) and 
which involve physical processes neglected in ASH such as 
ionization and radiative transfer. Solar values are used for 
the luminosity and rotation rate and the background stratifi- 
cation is based on a solar structure model. For further details 
on this simulation see Miesch et al. (2007). 

Near the top of the computational domain there is an in- 
terconnected network of downflow lanes which appears simi- 
lar to solar granulation but occupies a much larger scale (Fig. 
la). These giant cells are about 100 Mm across on average, 
compared to 30 Mm for supergranulation and 1 Mm for granu- 
lation. At high and mid latitudes, the downflow lanes exhibit 
intense cyclonic vorticity (counter-clockwise in the northern 




Figure 1: Flow features in a simulation of solar convection, (a) Radial velocity near the top of the computational domain 
(r = 0.98i?) with yellow tones denoting upflow and blue/black tones denoting downflow. This orthographic projection is 
tilted 35° toward the line of sight so the north pole is visible. (6) Differential rotation, (c) meridional circulation, and (d) 
specific entropy perturbations averaged over longitude and time (132 days). The meridional circulation is represented as 
streamlines of the mass flux and the entropy is normalized by the specific heat at constant pressure Cp = 3.5 x 10* erg g"-"^ 
K-i (after Miesch et al. 2007). 



hemisphere, clockwise in the southern), which is evident in 
animations of the convective patterns. 

The correlation time for the downflow network is 2-3 days, 
depending on latitude, but closer scrutiny reveals longer-lived 
structures (Miesch et al. 2007). At low latitudes, persistent 
downflow lanes occur which are oriented in a north-south di- 
rection and which persist for weeks to months. These are trav- 
eling convection modes which propagate prograde faster than 
the local rotation rate. Although the downflow network frag- 
ments deeper down, the north-south downflow lanes extend 
through most of the convection zone and play an important 
role in the maintcnence of mean flows (§3). 

The differential rotation achieved in this simulation is shown 
in Figure 16. It is roughly solar-like, with a monotonic de- 
crease in the angular velocity from equator to pole and rela- 
tively little variation with radius across the convection zone. 
However, the angular velocity contrast between the equator 
and ±60° is only about 50 nHz, compared to 90 nHz in the 
Sun. Other ASH simulations have achieved better agreement 
with helioseismic inversions, including a stronger angular ve- 
locity contrast and nearly radial angular velocity contours at 
mid latitudes (Miesch et al. 2006). 

Figure Ic shows the meridional circulation which is domi- 
nated by a single cell in each hemisphere, with poleward flow 
in the upper convection zone (amplitude 15-20 m s~^) and 
equatorward flow in the lower convection zone (amplitude 5- 
10 m s~^). This is roughly consistent with flux transport 
solar dynamo models and helioseismic determinations of the 
meridional circulation in the upper convection zone of the Sun 

(§!)■ 

Although the circulation pattern in Figure Ic is predom- 



inatly single-celled in each hemisphere, narrow counter-cells 
occur near the top and bottom boundaries. In light of our 
idealized boundary conditions (stress-free, impenetrable, im- 
posed specific entropy profile on the bottom and constant heat 
flux on the top) , the presence of these counter-cells should be 
interpreted with caution. ASH simulations which incorporate 
convective penetration into an underlying radiative zone gen- 
erally exhibit equatorward meridional circulation throughout 
the overshoot region (Miesch et al. 2000). Furthermore, an- 
gular momentum transport by supergranulation may induce 
a systematic poleward circulation in the surface layers via a 
process known as gyroscopic pumping (see eq. (1) below). 

Figure la, 6, and c illustrate convection, differential rota- 
tion, and meridional circulation. Before proceeding to §3, 
there is one more player in this drama to be introduced. 
Figure Id shows the mean specific entropy perturbation {S) 
relative to the spherically symmetric background stratifica- 
tion, averaged over longitude and time. The (S) profile ex- 
hibits prominent warm poles associated with thermal wind 
balance (§3). The corresponding temperature variation is ~ 
8K (Miesch et al. 2007, Fig. 6) which is more than five or- 
ders of magnitude smaller than the background temperature 
of 2.2 million K at the base of the convection zone. However, 
whereas the specific entropy perturbation increases monoton- 
ically with latitude throughout the convection zone, the tem- 
perature perturbation near the top of the shell typically peaks 
at the poles and the equator, with a minimum at mid-latitudes 
(Brun & Toomre 2002; Miesch et al. 2006). 
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Figure 2: Schematic diagram illustrating the complex nonlin- 
ear coupling between (a) convection, (b) differential rotation, 
(c) the mean specific entropy profile (S), and (d) the merid- 
ional circulation (see text). 

3 Maintenance of Mean Flows 

In order to understand how mean flows are maintained in our 
convection simulations and, by extension, in the Sun and in 
other stars, two equations are particularly enlightening. The 
first is derived from the zonal (longitudinal) component of the 
momentum equation, averaged over longitude and time: 

V- (p (v,„) C) = -V- {pr sin 9 (v'^K,)) (1) 

where C is the specific angular momentum 

C ^ r sin 9 {ilr sin 9 + (vcj,)) 

Equation ([1]) is expressed in a rotating spherical polar coordi- 
nate system with radius r, colatitude 9, and longitude and 
corresponding velocity components Vr,vg, and v^. The merid- 
ional velocity component is defined as = Vrr + ve9 where f 
and 6 are unit vectors in the r and 9 directions. Brackets <> 
denote averages over longitude and time, and primes indicate 
fluctuating variables for which the mean has been subtracted 
out, e.g. v'^ = — (vtj,). The rotation rate of the coordinate 
system is J7 and p is the (spherically symmetric) background 
density stratification. 

The second equation of particular importance is derived 
from the zonal component of the vorticity equation which is 
the curl of the momentum equation: 



9 d{S) 
2rCp 89 



(2) 



where (S) is the mean specific entropy perturbation as in Fig. 
Id, g is the gravitational acceleration, and Cp is the specific 
heat per unit mass at constant pressure. 

When deriving equations (1) and (2) we have assumed 
that the system in question (whether it be a simulation or a 
star) is in a statistically steady state. Furthermore, we have 
neglected viscous dissipation which is significant in some simu- 
lations but insignificant in stellar interiors where the Reynolds 
number (the ratio of inertial to viscous forces) is of order 10^^ 
or more. We have also neglected the Lorentz force which 
generally has a dissipative effect in the convection zone, sup- 
pressing rotational shear (Brun et al. 2004). In the tachocline, 
Lorentz forces associated with strong, localized toroidal field 
structures may divert or induce meridional and zonal flows. 



Equation (2) involves several additional assumptions, most 
notably that the Coriolis force operating on the differential 
rotation is large relative to the Reynolds stress. In other 
words, we have assumed that the Rossby number Ro ~ Lu/2n 
is much less than unity, where uj is the fluid vorticity in the 
rotating frame. The Rossby number is not specified in ASH 
simulations; rather, it is computed a posteriori from the sim- 
ulated flow field. Results generally indicate that the small 
Ro assumption implied by equation (2) is valid in the lower 
convection zone but breaks down near the surface where the 
amplitude of w peaks. 

When deriving equation (2) we have also assumed an ideal 
gas equation of state and that the background stratification 
is hydrostatic and adiabatic to lowest order. Both of these as- 
sumptions are well justified in the bulk of the solar convection 
zone but break down for r > 0.98i?. 

Subject to these reasonable assumptions, equation (1) states 
that the advection of angular momentum by the meridional 
circulation (vm) must balance the angular momentum trans- 
port due to the Reynolds stress (v'^V^. The Reynolds stress 

arises from correlations in the convective velocity components 
induced by the Coriolis force. Of particular relevance are the 
north-south oriented downflow lanes mentioned in §2; as hori- 
zontal flows near the top of the convection zone converge into 
these downflow lanes, eastward flows are diverted equator- 
ward by the Coriolis force and westward flows are diverted 
poleward. This produces an equatorward angular momentum 
transport which is mainly responsible for the prograde equa- 
torial rotation seen in Figure lb. 

Equation (2) states that gradients in the mean zonal ve- 
locity (as well as the mean angular velocity) parallel to the 
rotation axis are proportional to the latitudinal entropy gra- 
dient. This is the solar analogue of thermal wind balance 
which has been studied for decades within the context of geo- 
physics (e.g. Pcdlosky 1987). In the absence of latitudinal 
entropy gradients, equation (2) implies cylindrical rotation 
profiles in which angular velocity (and zonal velocity) con- 
tours are parallel to the rotation axis. This is a manifestation 
of the well-known Taylor-Proudman theorem. 

Insight into how the dynamical balances expressed by equa- 
tions (1) and (2) are achieved can be gained by considering 
Figure 2. Convection redistributes momentum and can es- 
tablish differential rotation (DR) and meridional circulation 
(MC) through the Reynolds stress (a — > 6, a — > d). It can also 
alter the mean entropy profile by means of the convective heat 
fiux which includes contributions from enthalpy fiux as well as 
kinetic energy fiux (o ^ c). Differential rotation can provide 
a negative feedback by shearing out convection cells (6 — > a). 
Meanwhile, modications of the specific entropy profile alter 
the buoyancy driving of the convective motions (c — > a) . The 
meridional circulation can in principle feed back directly on 
the convection as well {d a) but in practice this is negli- 
gible because the convective kinetic energy generally exceeds 
that contained in the MC by about two orders of magnitude. 

Differential rotation can induce meridional circulations 
through the Coriolis force (6 — d) while entropy variations 



can induce circulations through barochnicity {c ^ d), in par- 
ticular the d (S) 1 39 term on the right-hand-side of equation 
(2). Advection of angular momentum and entropy by the 
meridional circulation can in turn alter the (t;,^) and the {S) 
profiles {d ^ b, d ^ c). 

Equation (1) involves panels a, 6, and d of Figure 2, such 
that angular momentum transport by the convective Reynolds 
stress (a b) balances that by the meridional circulation 
{d ^ b). It is this balance that largely determines the merid- 
ional circulation pattern. As the Reynolds stress establishes 
a differential rotation, the meridional circulation adjusts such 
that equation (1) is satisfied. In previous, more laminar sim- 
ulations of convection, angular momentum transport by vis- 
cous diffusion upset the balance expressed by equation (1) and 
circulation patterns were qualitatively different from that in 
Figure Ic, multi-celled in latitude and radius (Miesch et al. 
2000, 2006; Brun & Tooomre 2002). 

This simplified picture is complicated by the nonlinear 
feedback mechanisms illustrated in Figure 2 and the addi- 
tional constraint provided by thermal wind balance, equation 
(2). Equation (2) only explicitly involves the differential ro- 
tation and specific entropy profile (6, c), but convection and 
meridional circulation (a, d) both play an essential role. In 
order to demonstrate this, we consider now two illustrative 
scenarios by which thermal wind balance may be achieved. 
We'll refer to these as mechanical and thermal forcing (Mi- 
esch et al. 2006). Both cases begin with convection (a) but 
they each circle the diagram in Figure 2 in an opposite sense, 
clockwise and counter-clockwise respectively. 

In the mechanical forcing scenario, the convective Reynolds 
stress establishes a differential rotation (a ^ b) which in- 
duces a meridional circulation via the Coriolis force (6 — > c). 
The advection of entropy by these circulations then alters the 
background stratification {S) (c ^ d). In the subadiabatic 
portion of the tachocline, this would tend to establish warm 
poles and the Coriolis-induced circulations would cease when 
thermal wind balance, equation (2), is achieved (Rcmpel 2005; 
Miesch et al 2006). However, in the convection zone the sense 
of the induced circulation is such that the resulting latitudinal 
entropy gradient would be equatorward rather than poleward 
(cool poles). Thermal wind balance would be achieved even- 
tually, but the rotation profile woidd not be solar-like. 

Now consider thermal forcing in which a latitude-dependent 
convective heat flux establishes warm poles (o — > c) which 
in turn induce a meridional circulation through barochnic- 
ity (c — *■ d). Redistribution of angular momentum by this 
baroclinic circulation then establishes a differential rotation 
{d ^ b). This will 'work' in the sense that the induced DR 
will evolve until equation (2) is satisfied. However, in the ab- 
sence of the Reynolds stress, the meridional circulation would 
tend to conserve angular momentum, accelerating the poles 
relative to the equator. This follows from equation (1); if the 
right-hand-side is zero and V- (p {vm)) = as required by the 
anelastic approximation, then (v^) - VC = 0, so £ would be 
constant along circulation streamlines. A solar-like differen- 
tial rotation profile can only be achieved if the Reynolds stress 
also contributes, transporting angular momentum toward the 



equator (a b). 

Thus, the convective Reynolds stress is needed to produce 
a prograde equatorial rotation relative to higher latitudes as 
in the Sun while barochnicity is needed to break the Taylor- 
Proudman preference for cylindrical rotation profiles. Such 
barochnicity arises from a latitude-dependent convective heat 
flux but thermal coupling to the tachocline may also play an 
important role (Rempel 2005; Miesch et al. 2006). Much has 
been learned but much remains to be done and convection 
simulations will continue to provide valuable insight into how 
differential rotation and meridional circulation are established 
and maintained in the Sun and in stars across the HR dia- 
gram. 
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DeRosa, Matthew Browning, Benjamin Brown, Nicholas Feath- 
erstonc, Kyle Augustson, and Nicholas Nelson. Funding was 
provided by NASA through grant NNG05G124G of the He- 
liophysics Theory Program. The simulations were carried 
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